Detecting the breached pair phase in a polarized ultracold Fermi gas 
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We propose a method for the experimental detection of a new quantum phase, the breached pair 
state, in a strongly interacting ultracold Fermi gas with population imbalance. We show that through 
the time-of-flight Raman imaging, the presence of such a phase can be unambiguously determined 
with a measurement of the momentum-space phase separation of the minority spin component. To 
guide the experimental efforts, the momentum-space density profiles are calculated under typical 
experimental conditions. 



The search for new quantum phases has long been a 
persistent goal of many-body physics. Ultracold atomic 
gas, with its remarkable controllability, has provided 
a powerful platform to realize and probe such phases. 
Among the various many-body phases, a few have been 
realized in this system, including, for instance, the Bose 
or the Fermi superfluid states 0| , and the Mott insulator 
state in an optical lattice 0|- Here, we propose a method 
for the experimental detection of a new quantum phase, 
the breached pair state, in a strongly interacting gas of 
fermionic atoms with population imbalance. Such a sys- 
tem, with its recent experimental realization [a, kl la : has 
raised strong interest B H H E3 El Hill H 13 
Hfil ll^ l . In particular, a quantum phase transition from 
a BCS superfluid to a normal state has been observed 



The breached pair (BP) state, proposed in Refs. 
H @ and studied in detail in [ill [H [H Q, is a 
non-BCS superfluid phase with gapless fermionic excita- 
tions. A unique feature of this state is its phase separa- 
tion in the momentum space, although in real space it is 
a homogeneous, polarized superfluid. This feature also 
distinguishes the breached pair phase from a simple mix- 
ture of molecular condensate and fermionic atoms, where 
there is no phase separation in the momentum space. As 
illustrated in Fig. 1, there are two kinds of BP states, 
with one or two Fermi surfaces for the excess fermions 
0E1j respectively. To distinguish the two scenarios, we 
simply label them the BPl and the BP2 states. Except 
for the number of Fermi surface, the properties of these 
two states are very similar. The quasi-particle excitations 
in both of the BP phases have non-monotonic dispersion 
relations 0, 0, E3 • There exists a well-defined quantum 
phase transition from the BCS superfluid to either one 
of the two phases. Such a phase transition is not charac- 
terized by symmetry breaking, but rather by the change 
in the topology of the Fermi surface 0, Q] . 

The BP phases have been predicted to exist in the- 
ory, but it is hard to detect them experimentally. These 
phases show polarized superfluidity. However, even if one 
can confirm polarized superfluidity in the future experi- 
ments, it is not an unambiguous signal for the BP phase. 
Many phases can have polarized superfluidity. In par- 
ticular, as we have shown in the conventional BCS 
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FIG. 1: (a) The level configuration for the 6 Li atoms near the 
Feshbach resonance point Bo = 834G. An impulsive n Raman 
pulse can selectively transfer the atoms from the lowest level 
(II)) to the highest level (|6)) in the ground-state manifold, 
(b) Schematic momentum-space configurations of fermions in 
the BPl (top) and the BP2 (bottom) phases. The shaded 
regions with dots are filled with the unpaired fermions of the 
majority spin component, while the regions labelled 'SF' are 
filled with the gapless superfluid. (c) Schematic momentum- 
space density profiles of the minority fermions in the BPl 
(left) and the BP2 (right) phase. 



state in the polarized gas also gives polarized superflu- 
idity at any finite temperature (which is necessarily the 
case for experiments) due to the polarized quasi-particle 
excitations. So, it is hard to distinguish the more exotic 
phases from the conventional thermal BCS superfluid in 
experiments. 

To overcome this difficulty, we propose a method to de- 
tect the BP phases by measuring the momentum-space 
density profile of the minority spin component. This 
method has the following desirable features: first, it can 
give an unambiguous signal for the BP phases. From 
the basic physical picture of the BP states shown in Fig. 
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1(b), we see already that there exists a finite region in 
the momentum-space with no minority atoms. This re- 
gion, enclosed by the Fermi surface(s), is occupied only 
by the majority spin component. Therefore, if one mea- 
sures the momentum-space density profiles, the profile 
for the minority atoms should be non-monotonic with a 
dip in the momentum distribution. This dip, if observed, 
is an unambiguous signature of the BP phase as no other 
phases have that character. Second, this method can also 
unambiguously distinguish different kinds of BP states. 
In Fig. 1(c), we show schematically the characteristic 
momentum distributions for the minority atoms associ- 
ated with different BP phases. The different structures 
of the dips in the momentum distributions clearly tell 
the nature of the corresponding phases. More detailed 
calculations of the momentum space profiles later in the 
paper will further confirm the feasibility of this detection 
method. 

For non-interacting atomic gases, the time-of-flight 
imaging provides a way to measure the momentum-space 
density profile. The absorption image of the expanded 
cloud is proportional to the initial atomic momentum dis- 
tribution (after column integration). However, to realize 
the BP states, the atoms need to be strongly interacting 
in the trap. One therefore has to devise a method to 
suddenly turn off this strong interaction right before the 
expansion. For 40 K atoms, such an interaction turn-off 
can be achieved with a fast switch of the magnetic field, 
as has been demonstrated in the recent experiment |l8| . 
The experiments on the polarized Fermi gases are done so 
far with 6 Li atoms, whose inter-atomic interaction, due to 
the large resonance width of 6 Li at Bq = 834G, is hard to 
be turned off through the magnetic filed switching with- 
out introducing significant influence on the momentum 
distribution of the atoms. 

To overcome this problem, we propose a method to 
tune the atomic interaction almost instantaneously by 
applying an impulsive Raman pulse consisting of two co- 
propagating laser beams. For 6 Li atoms close to a Fes- 
hbach resonance (the magnetic field ranging from 600 G 
to 1000 G), the energy eigen-levels in the ground-state 
manifold are shown in Fig. 1(a). The three lower (or up- 
per) levels differ mainly by their respective nuclear spins, 
while the states between the two sets have different elec- 
tronic spins. The Feshbach resonance at Bo = 834G res- 
onantly couples the lowest two energy levels (|J.) and ||)) 
in the ground state manyfold. Near this resonance point, 
the intra-set level splitting (e.g., between the states |J.) 
and |t)) is about 70 MHz, while the inter-set splitting 
(e.g., between the states ||) and |6)) is around 2.5 GHz. 
We assume that the atoms are initially in the levels |t) 
and |t), with spin J, as the minority component. Right be- 
fore the time-of-flight imaging, we apply an impulsive Ra- 
man pulse to fully transfer the minority component from 
the state |t) to |6) (see Fig. 1(a)). As the level |6) is de- 
tuned from |t) by 2.5 GHz, which is larger than the width 



of the Feshbach resonance at B (2^bW//i ~ 0.84 GHz, 
where W = 300G is the width of the Feshbach resonance 
at Bo = 834G), it does not participate in the resonant 
coupling and thus only weakly interacts with the spin |t) 
component. Since the gas is dilute and the atoms in the 
two different states |t) and |6) can be quickly separated, 
as a good approximation, the momentum distribution of 
the 1 6) component has a negligible variation during the 
free expansion and can be recorded with the absorption 
image. As the Raman pulse consists of co-propagating 
laser beams, the momentum distribution of the atoms 
are not changed during the transfer. Therefore, using 
the method outlined above, one effectively measures the 
in-trap momentum-space density profile for the minority 
atoms. 

We add a few remarks on this detection method be- 
fore going on to the theoretical calculation. Firstly, in 
this scheme the Raman pulse is impulsive with its dura- 
tion significantly smaller than h/ /j,, where fj, is the char- 
acteristic chemical potential of the gas. The dynamics 
of the system during such a short interval is negligible, 
and thus the Raman transfer faithfully keeps track of the 
momentum distribution of the atoms. As a result of the 
impulsiveness of the Raman pulse, this method does not 
resolve the energy spread in the atomic gas, and is there- 
fore different from the Bragg or the r.f. spectroscopy 
[lH Efij l . The two methods are actually complementary 
to each other, as the spectroscopy measures the energy 
spread while the Raman imaging detects the momentum 
distribution. Secondly, an impulsive Raman pulse with- 
out resolving the energy spread may have a practical ad- 
vantage as it is insensitive to small variation of the en- 
ergy splitting between the levels |t) and |6) caused by the 
fluctuation of the magnetic field. Finally, we can also let 
the two laser beams of the Raman pulse propagate along 
different directions (say, along x and — x axis). Such a 
pulse will not change the momentum distribution in the 
y — z plane, while imprinting a momentum kick along the 
x direction, which helps to quickly separate the atomic 
components |t) and |6), and thus further reduces the in- 
fluence of interaction during the expansion. 

To guide the experimental effort, in the following we in- 
clude some calculations on the momentum-space density 
profile of the minority atoms. The calculations show that 
on the BEC side of the resonance |2l| , a non-monotonic 
momentum profile with a dip in the center should be vis- 
ible, even if we take average in real space over all the 
phases in the trap [22T | and perform a column integration 
in the momentum space. Our calculation method |l3j | is a 
gen eralization of the self-consistent GoG diagram scheme 
[23j to the case with unequal spin population, which 
at zero temperature reduces to the mean-field approach 
[llHl2ill3Lll4l | (the latter recovers the crossover theory in 
the equal population case ). We use the local density 
approximation for a three-dimensional (anisotropic) har- 
monic trap, and pick up the phase corresponding to the 
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global minimum of the thermodynamic potential. The 
detailed formalism has been described in Refs. 0|. The 
properties of the system are universally characterized by 
three dimensionless parameters, the population imbal- 
ance (3 = (iVf — iV;) / (iVj + iVjJ, the dimensionless in- 
teraction strength kpa s , and the temperature T/Tp , 
where kp = y // 2mEp/h 2 and Tp = Ep/ks are defined 
through the Fermi energy Ef at the center of the trap 
for a non-interacting Fermi gas with equal spin popula- 
tions under the local density approximation. ^From this 
definition, Ep = (3Nh 3 uj x ujyUj z )3 , where N = + N± is 
the total particle number, (i = x, y, z) is the trapping 
frequency along the ith direction. The atomic scattering 
length a s is connected with the magnetic field detuning 
AB = B — Bo through a s = ab g (l — W/AB), where Ob s 
is the background scattering length and W is the width 
of the Feshbach resonance. 

The phase diagram of this system has been shown in 
fl3|. according to which the BP1 phase appears on the 
BEC side of the resonance with (kpas)^ 1 > 0.5 at low 
temperature. To calculate the momentum space den- 
sity profile, however, we need to average over all the 
phases in the trap |24|. Under the local density approx- 
imation, the atom number distribution n\^ r depends on 
both the atomic momentum k and the position in the 
trap r. The momentum-space density profile is given 
by ?i k = / n kr d 3 r. At (fc F a s ) _1 = 2, j3 = 0.8, and 
at zero temperature, this density profile rtk is plotted 
in Fig. 2 for the minority spin component (spin j), to- 
gether with its column integration n'^ = J n^dkz . It is 
clear that a non-monotonic behavior should be visible 
for this density profile, even after the column integration 
[25| . The dip is smaller than the one in the schematic 
figure 1(c) as one needs to average over different phases 
here in the trap. With the above parameters, there is a 
large region of the BP1 phase near the center of the trap 
(there is still a small BCS superfluid core right at the cen- 
ter), which is surrounded by a normal gas of the major- 
ity spin component. To achieve the interaction strength 
{kpa s )~ l = 2, one can take a moderately negative de- 
tuning with AB — —80 G for the 6 Li atoms, where the 
system still has a long lifetime. This parameter then cor- 
responds to N = 1.5 x 10 5 atoms in a trap with typical 
trap frequencies ui x = oj y ~ 110Hz and to z ~ 2ZHz 3] 
(or with N = 1.2 x 10 6 atoms in a slightly weaker trap 
with the trap frequencies Ui reduced by a factor of 2). 

To check stability of the signal for the momentum- 
space phase separation even at finite temperature, and to 
verify that this signal is indeed uniquely connected with 
the BP phase, we calculate the momentum-space den- 
sity profiles nk at different interaction strengths (kpa s ) 
and temperatures (T/Tp), and with different popula- 
tion imbalance parameters {(5). The results are shown 
in Fig. 3, from which several observations can be made. 
First of all, this signal (the non-monotonic profile nk) 
remains stable at finite but low temperatures. For in- 




FIG. 2: (a) The momentum-space density profile (only 
dependent on k — |k|) and (b) its column integration n' k 
for the minority spin component. The unit of the density is 
riF = k%/(3ir 2 ), where &f has been specified in the text. 
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FIG. 3: The momentum-space density profile nk at different 
temperatures, interaction strengths, and population imbal- 
ance. The solid lines represent the zero-temperature distri- 
butions, and the dashed ones represent the distributions at 
T = 0.2Tf. The interaction strength (kpa 3 ) and the popu- 
lation imbalance /3 are specified in the figure. Note that the 
hump at small momenta for the solid line in panel (d) indi- 
cates the existence of a normal mixed phase (NM) in the trap 







stance, at T = Q.2Tp (which corresponds a real temper- 
ature T ~ 50ni^ for N = 1.5 x 10 5 atoms in a trap with 
lu x = uj y ~ 110-ffz and oj z ~ 23Hz), the dip in the pro- 
file only becomes a bit shallower. At significantly higher 
temperatures, the BP state is eventually replaced by a 
normal phase, and the dip signal disappears. Secondly, 
when one goes slightly further to the BEC side [26] , such 
as with (kpas)^ 1 — 3, the dip signal becomes more pro- 
nounced, which corresponds to a larger BP phase region 
at the trap center. On resonance and on the BCS side, 
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the dip disappears (the dip first disappears on the BEC 
side where the BP phase ceases to exist in the trap [l3|h 
Finally, as one varies the interaction strength kpa s or the 
population imbalance /?, the dip size varies, which cor- 
relates with the variations in the size of the BP region 
in the phase diagram. This shows that the dip signal is 
indeed uniquely connected with the BP phase. 

Before ending the paper, a few remarks are in order. 
Firstly, in our calculation, only the BP1 phase can be 
stabilized in this system, and it only appears on the 
BEC side of the resonance |2j|. Our proposed detection 
method, however, can detect both the BP1 and the BP2 
phases, and an experimental detection surely can scan 
the whole region across the Feshbach resonance. Similar 
to the real-space density profile [3J, |^, |5j , a momentum- 
space density profile will help us to understand a lot 
about this strongly interacting system with various com- 
peting phases, whether the experiment confirms the the- 
oretical prediction here or finds new surprise. Secondly, 
we use a single impulsive Raman pulse here to fast tune 
the atomic interaction. If one applies two consecutive 
impulsive Raman pulses propagating along different di- 
rections, as shown in Ref. |28j . one can reconstruct the 
full real-space or momentum-space correlation function 
through Fourier sampling with the laser phase gradient. 
The full correlation function can reveal very detailed in- 
formation of the system, e.g., the tiny spatial structures 
in the system. With such a capability, one may find some 
other interesting quantum phases, such as the FFLO 
state [la,|22ll2fj, if the latter does show up in this atomic 
gas. 
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